Search for: All records

We introduce an individually fitted screened-exchange interaction for the time-dependent Hartree–Fock (TDHF) method and show that it resolves the missing binding energies in polymethine organic dye molecules compared to time-dependent density functional theory (TDDFT). The interaction kernel, which can be thought of as a dielectric function, is generated by stochastic fitting to the screened-Coulomb interaction of many-body perturbation theory (MBPT), specific to each system. We test our method on the flavylium and indocyanine green dye families with a modifiable length of the polymethine bridge, leading to excitations ranging from visible to short-wave infrared. Our approach validates earlier observations on the importance of inclusion of medium range exchange for the exciton binding energy. Our resulting method, TDHF@vW, also achieves a mean absolute error on a par with MBPT at a computational cost on a par with local-functional TDDFT. 

We develop an efficient approach to evaluate range-separated exact exchange for grid- or plane-wave-based representations within the generalized Kohn–Sham–density functional theory (GKS–DFT) framework. The Coulomb kernel is fragmented in reciprocal space, and we employ a mixed deterministic-stochastic representation, retaining long-wavelength (low-k) contributions deterministically and using a sparse (“fragmented”) stochastic basis for the high-k part. Coupled with a projection of the Hamiltonian onto a subspace of valence and conduction states from a prior local-DFT calculation, this method allows for the calculation of the long-range exchange of large molecular systems with hundreds and potentially thousands of coupled valence states delocalized over millions of grid points. We find that even a small number of valence and conduction states is sufficient for converging the HOMO and LUMO energies of the GKS–DFT. Excellent tuning of long-range separated hybrids (RSH) is easily obtained in the method for very large systems, as exemplified here for the chlorophyll hexamer of Photosystem II with 1320 electrons. 

We present a highly efficient method for the extraction of optical properties of very large molecules via the Bethe–Salpeter equation. The crutch of this approach is the calculation of the action of the effective Coulombic interaction, W, through a stochastic time-dependent Hartree propagation, which uses only ten stochastic orbitals rather than propagating the full sea of occupied states. This leads to a scaling that is at most cubic in system size with trivial parallelization of the calculation. We apply this new method to calculate the spectra and electronic density of the dominant excitons of a carbon-nanohoop bound fullerene system with 520 electrons using less than 4000 core hours. 

Excitonic chromophore aggregates have wide-ranging applicability in fields such as imaging and energy harvesting; however their rational design requires adapting principles of self-assembly to the requirements of excited state coupling. 

Molecular aggregates with long-range excitonic couplings have drastically different photophysical properties compared to their monomer counterparts. From Kasha's model for one-dimensional systems, positive or negative excitonic couplings lead to blue or red-shifted optical spectra with respect to the monomers, labeled H-and J-aggregates, respectively. The overall excitonic couplings in higher dimensional systems are much more complicated and cannot be simply classified from their spectral shifts alone. Here, we provide a unified classification for extended 2D aggregates using temperature dependent peak shifts, thermal broadening, and quantum yields. We discuss the examples of six 2D aggregates with J-like absorption spectra but quite drastic changes in quantum yields and superradiance. We find the origin of the differences is, in fact, a different excitonic band structure where the bright state is lower energy than the monomer but still away from the band edge. We call this an “I-aggregate.” Our results provide a description of the complex excitonic behaviors that cannot be explained solely on Kasha's model. Furthermore, such properties can be tuned with the packing geometries within the aggregates providing supramolecular pathways for controlling them. This will allow for precise optimizations of aggregate properties in their applications across the areas of optoelectronics, photonics, excitonic energy transfer, and shortwave infrared technologies.